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METHOD AND APPARATUS FOR CALIBRATING BASE STATION LOCATIONS 
AND PERCEIVED TIME BIAS OFFSETS IN AN ASSISTED GPS TRANSCEIVER 



Technical Field 

5 The present invention relates generally to a satellite positioning receiver, and 

more particularly to a Global Positioning System (GPS) receiver in an assisted GPS 
handset operating in conjunction with one or more cellular base stations (BSs). 
Background Art 

The Global Positioning System (GPS) is a satellite-based positioning system 

10 developed by the U.S. Department of Defense to give accurate positional information to 
a GPS receiver anywhere in the world. A properly equipped GPS receiver may therefore 
be used in any setting in which a positional fix is desired, and typically yields positional 
coordinates in three dimensions. 

GPS is enabled by a satellite orbital constellation made up of twenty-four or more 

15 satellites orbiting the earth in twelve-hour orbits. The satellites are arranged in six orbital 
planes, each containing four satellites. The orbital planes are spaced sixty degrees apart 
and are inclined approximately fifty-five degrees from the equatorial plane. This 
constellation ensures that from four to twelve satellites will be visible at any time at any 
point on earth with a clear view of the sky. 

20 A three-dimensional position fix may be determined if signals are being received 

from four or more GPS satellites. The received satellite signals each contain an identifier 
unique to that particular satellite. These identifier codes are commonly called Gold codes 
or pseudorandom noise (PRN) codes and allow a GPS receiver to discriminate between 
signals from different satellites. Also contained within the signals are satellite ephemeris 

25 data containing information such as an orbital configuration and a satellite time. (All 
GPS satellite signals contain a common, synchronized GPS time.) This satellite time 
signal allows a GPS receiver to detect a time of receipt and therefore measure a transit 
time of the signal. In turn, the transit time enables a GPS receiver to determine a distance 
(termed a pseudorange) to the satellite. Four GPS pseudoranges allow a location on the 

30 Earth's surface to be determined. 

There are many uses for portable, hand-held GPS receivers including use in an 
assisted GPS handset. Because of processing time constraints and incomplete satellite 

i 
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coverage, handset-based GPS receivers may operate in conjunction with one or more 
BSs. As described in related art, the base station (BS) can be employed to relieve GPS 
receivers from certain processing duties, improve the accuracy of the resulting location 
solution, and improve fix coverage in areas that are marginal for GPS signal tracking 
5 (e.g., in an urban canyon or indoors). The serving BS (i.e., the BS that is responsible for 
maintaining communication with the handset) can provide the mobile handset with 
satellite Doppler information (U.S. Patent No. 4,445,118 to Taylor, et al.) and/or satellite 
ephemeris data (U.S. Patent No. 5,365,450 to Schuchman, et al.) to reduce the time 
required to acquire the GPS signals and derive a location fix within the handset. In 

10 addition, differential corrections can be broadcast to GPS receivers in handsets. These 
corrections, if received with minimum latency, can remove most of the error associated 
with each pseudorange measurement. Finally, additional range measurements can be 
derived from the arrival times associated with the transmissions of pilot signals from 
each BS, which can be used to augment the GPS-only solution and provide position fixes 

15 in situations where GPS cannot. One type of information that is commonly relayed in the 
related art is locations of nearby BSs. BS locations may be useful to the mobile GPS 
receivers to give a starting position for a positional fix and may reduce the signal search 
space and therefore reduce the acquisition time for the GPS signals. A second type of 
information that is commonly relayed is a time bias offset associated with each nearby 

20 BS. This offset is a bias between the time maintained by the BS and GPS time. 

FIG. 1 shows a related art mobile GPS receiver in a cellular handset operating in 
conjunction with a single BS. Each related art cellular handset receives a transmission 
from the BS containing the BS position embedded therein along with a time bias offset b. 
The BS position may be expressed in Cartesian coordinates as y BS . A drawback of the 

25 related art approach is that related art handset-based GPS receivers cannot make use of 
TOA measurements from the network infrastructure and calculate a positional fix based 
solely on TOA measurements until the required data (x^ y^, b) is received from the BS. 

What is needed, therefore, is a method for determining BS locations and time bias 
offsets in a handset-based GPS receiver without necessarily having to receive the BS 

30 locations and time bias offsets from the infrastructure. 



2 
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Brief Description of Drawings 

FIG. 1 shows related art mobile GPS receivers operating in conjunction with a 
base station; 

FIG. 2 illustrates a distance determination utilizing Cartesian coordinates; 
5 FIG. 3 shows a typical mobile station to base station geometry illustrating that the 

mobile station is moving relative to the base station; 

FIG. 4 shows a flowchart of a first embodiment of a method for calibrating base 
station locations and time bias offsets; 

FIG. 5 shows a flowchart of a second embodiment of a method for calibrating base 
10 station locations and time bias offsets; and 

FIG. 6 shows a GPS receiver that can be used for calibrating base station locations 

and time bias offsets. 



Disclosure of the Invention 

A method for calibrating a base station (BS) location and/ or perceived time bias 

15 offset in a mobile station (MS) with a satellite positioning receiver, such as a GPS 

receiver, is provided according to a first embodiment. The method includes the steps of 
collecting a number of TOA measurements arriving at the MS from a BS, simultaneously 
collecting a corresponding number of GPS position measurements of a position of the 
GPS receiver, and computing the BS location and/ or perceived time bias offset using the 

20 collected TOA and GPS position measurements. 

A method for calibrating a BS location and/or perceived time bias offset in an MS 
with a satellite positioning receiver, such as a GPS receiver, is provided according to a 
second embodiment. The method includes the steps of collecting at least three TOA 
measurements TOA,, TOA 2 and TOA 3 arriving at the MS from a BS, wherein a TOA 

25 measurement is capable of being expressed as TOA, = b' + V[(x-x,) 2 + (y-y,) 2 ], with 
variable b' being a perceived time bias offset as seen by the MS and (x, y) being 
unknown BS coordinates. The method simultaneously collects at least three GPS position 
measurements of a position of the GPS receiver, wherein a position measurement is 
denoted as x, and y, as a function of time. The method then subtracts a first TOA 

30 measurement from two subsequent TOA measurements to create two TOA 

measurement differences TOA 31 and TOA^ and substitutes a distance variable d, for the 
unknown BS coordinates x, y in the two TOA measurement differences TOA 31 and TOA^ 
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to produce a quadratic equation for the distance variable d u where the distance variable 
is given by d 1 = V[(x-x,) 2 + (y-y t ) 2 ]- The method then solves the quadratic equation for the 
distance variable d t to obtain an estimated range d^ solves for the unknown BS 
coordinates x, y using the estimated range d e , and obtains the perceived time bias offset 
5 b' by subtracting the estimated range d e from TOA ir TOA,, or TOA 3 . 

An apparatus for calibrating a BS location is provided according to a third 
embodiment. The apparatus includes an MS with a satellite positioning receiver, a TOA 
measurement storage capable of storing a plurality of TOA measurements, a GPS 
receiver position fix storage capable of storing a plurality of position fixes, a BS 

10 coordinates solution storage, a time bias offset storage, and a processor. The processor 
communicates with the satellite positioning receiver to receive signals from positioning 
satellites and the BS and further communicates with the TOA measurement storage, the 
satellite positioning receiver position fix storage, the BS coordinates solution storage, and 
the time bias offset storage. The processor is capable of collecting and storing TOA 

15 measurements, collecting and storing a corresponding number of satellite position 
measurements, subtracting a first TOA measurement from subsequent TOA 
measurements to create TOA differences, substituting a distance variable d, for the 
unknown BS coordinates x, y in the TOA differences to produce a quadratic equation, 
solving the quadratic equation for the distance variable d, to obtain an estimated range 

20 d e , solving for the unknown BS coordinates x, y, obtaining the perceived time bias offset 
by subtracting the estimated range d e from the TOA measurements, linearizing TOA 
equations to produce linearized TOA equations, generating a vector of measurement 
residuals TOA res , computing a weighted least squares solution x^, of the TOA equations, 
using the weighted least squares solution x so , to assess an accuracy of the linearized TOA 

25 equations, and stopping the linearizing when solution error variances of the weighted 
least squares solution x., become less than a solution error variance threshold. 

Once the BS calibration is essentially complete, the perceived time bias offset b' as 
seen by the MS is periodically re-determined and associated with each MS location or BS 
service area. Over time, this association permits construction of a temporal map, or 

30 database, of perceived time bias offsets versus position over the area served by a BS. The 
perceived time bias offset b' represents the time bias offset b as affected by spatial error 
from multipath in addition to the fixed error induced by cable delays, internal BS 
harnessing, software processing, clock source distribution, and other factors. This map 
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can be derived and maintained for many BSs by the cellular infrastructure and broadcast 
to other MSs (whether or not GPS-equipped) to improve their location accuracy. 

FIG. 2 illustrates a distance determination utilizing Cartesian coordinates. 
Although Cartesian coordinates are used here, other coordinates (e.g. polar) can be used 
5 to determine distance. In the figure, the distance d> from the MS to the BS can be 
determined by the formula d, = V[(x, - x 0 ) 2 + (y, - y 0 ) 2 ]. This formula assumes that the 
earth is locally flat over the distance between the BS and MS, an excellent approximation 
over the range of distances expected (i.e., less than twenty km). An MS as used herein 
can be a cellular handset, a vehicle-based communication transceiver, or any other type 

10 of mobile radio transceiver. 

FIG. 3 shows a typical MS to BS geometry having an MS with a satellite 
positioning receiver that is moving relative to the BS. Because of the movement of the 
MS, it is at position MS, at time t = 1, at position MS 2 at time t = 2, at position MS 3 at time 
t = 3, etc. Therefore, the distance to the BS may change, as shown by distances d u d 2 , and 

15 d 3 . 

FIG. 4 shows a flowchart 400 of a first embodiment of a method for calibrating BS 
locations and perceived time bias offsets. In step 402, MS TOA measurements are 
collected. The derivation of the BS coordinates x, y and the perceived time bias offset 
b' may be done through use of at least three received TOA measurements from the BS. 
20 In use, the MS may receive multiple TOA measurements from the BS. A first 

derived TOA measurement is: 

TOA, = b' + V[(x - x t ) 2 + (y - y,) 2 ] (1) 
where TOA t is a TOA measurement derived by the BS at a time t = 1, where x t and y t are 
the known GPS position fix coordinates of a mobile GPS receiver at a time t (i.e., x 2 and y t 
25 are x t and y t at the time t = 1), where b' is the unknown perceived time bias offset of the 
BS, and where (x, y) are the unknown BS coordinates. Equation (1) therefore represents a 
single equation having three unknowns. 

Equation (1) is linear in the perceived time bias offset b' variable, but nonlinear in 
the BS coordinates x, y. Note that the TOA measurement TOA t and the mobile GPS 
30 receiver coordinates x t , y t are time variant. 

At a later time t = 2, denoted by the subscript 2, a second TOA measurement is 
derived from the BS: 

TOA 2 = b' + V[(x - x/ + (y - y 2 ) 2 ] (2) 
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Consideration of equations (1) and (2) provides two equations with three unknowns. 
Derivation of a third TOA is required to generate estimates of x, y, and b'. At yet a later 
time t = 3, denoted by the subscript 3, a third TOA measurement is derived from the BS: 

TOA 3 = b' + V[(x - x 3 ) 2 + (y - y 3 ) 2 ] (3) 

5 Note that a single perceived time bias offset term b' is used in equations (l)-(3). The 

change in b' due to movement of the GPS receiver is assumed sufficiently small so that b' 
can be assumed a constant. Once the BS location has been accurately determined and 
verified, this assumption can be relaxed, and the perceived time bias offset b' can be 
determined periodically. 

10 In step 406, a corresponding number of GPS receiver position measurements are 

generated by the GPS receiver. These take the form of x t and y t , where t is the time of 
measurement. Therefore, as the GPS receiver moves, it generates distinct positional fixes 
(x^y,), (x 2 ,y 2 ), and (x 3 ,y 3 ), as shown in FIG. 3 and equations (l)-(3). 

In step 408, the singular geometry of the TOA measurement set is tested. This is 

15 done through use of a matrix M, given as: 

I (Xa-x,) (y 3 -yi) I 

I (X.-X,) (y 2 -yi) I 

The determinant of the matrix M above (which must be nonzero to solve the equations) 

20 is proportional to the area enclosed by the three MS positions. 

More generally, any near singularity condition that occurs when inverting the 
two-dimensional matrix Mis an indication that the measurement set is not sufficient to 
determine the position of the BS, however coarsely. If the singularity of the received 
TOA measurements is inadequate, the method may acquire more TOA measurements. 

25 In step 410, the TOA meastirements and the GPS receiver position measurements 

are use to compute the (unknown) BS location coordinates x, y. Note that a significant 
geometry change must be induced by the motion of the MS before the equations can be 
solved. For example, if the MS was stationary, the two resulting linear equations in x and 
y would be linearly dependent and so could not be solved. 

30 Since computation of the BS coordinates is more difficult than the computation of 

the perceived time bias offset b' and involves nonlinear processing, it is addressed first. 
It should be understood that alternately the perceived time bias offset b' could be 

6 
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derived first and the BS location could then be derived knowing the perceived time bias 
offset. 

Determination of an initial solution may be performed as follows. Equation (1) 
may be subtracted from both equations (2) and (3), producing the TOA differences: 
5 TOA 21 = TOA 2 - TOA, = V[(x - + (y - y 2 ) 2 ] - V[(x - x,) 2 + (y - y,) 2 ] (4) 

TOA 31 = TOA3 - TOA, = V[(x - X3) 2 + (y - y 3 ) 2 ] - V[(x - x,) 2 + (y - y,) 2 ] (5) 
where the perceived time bias offset b' drops out due to the subtraction. Now, the 
distance variable d, from the MS to the BS at time 1 (see FIG. 2) may be defined as: 

d 1 = V[(x-x 1 ) 2 + (y-y 1 ) 2 ] (6) 
10 The -V[(x - x,) 2 + (y - y,) 2 ] terms at the right of equations (4) and (5) are equal to -d, by 
definition. Therefore, by adding +d, to both sides of equations (4) and (5), they may be 
rewritten as: 

TOA 31 + d, = V[(x - x 3 ) 2 + (y - y3 ) 2 ] (7) 
TOA 21 + d 2 = V[(x - x 2 ) 2 + (y - y 2 ) 2 ] (8) 

1 

15 Squaring equations (7) and (S) results in: 

TOA 31 2 + 2d,TOA 31 + d, 2 = (x - x 3 ) 2 + (y - y 3 ) 2 (9) 
TOA 2I 2 + 2d,TOA 2 , + d, 2 = (x - x,) 2 + (y - y 2 ) 2 (10) 

Subtracting the term d, 2 from both sides of equations (9) and (10) gives the desired 

results: 

20 TOA 3I 2 + 2d 1 TOA 31 = -2(x 3 - x,)x - 2(y 3 - y ,)y + k 3 - k, (1 1) 

TOA,, 2 + 2d l TOA 21 = -2(x 2 - x,)x - 2(y 2 - y,)y + k, - k, (12) 

where 

K = x, 2 + y, 2 
k2 = x 2 2 + y 2 2 

25 i<3 = x 3 2 + y 3 2 

Equations (11) and (12) represent two linear equations, which can be solved for x 
and y in terms of the unknown distance variable d, (the distance from the MS to the BS at 
time t = 1), provided that the matrix M given previously is invertable. If the matrix M 
cannot be inverted, additional MS positions and TOA measurements must be collected 
30 before the location can be accurately determined. Given x and y in terms of the distance 
variable d, (i.e., equation (6)), the distance variable d, can be substituted for the term V[(x 
_ x ) 2 + (y . yj ) 2 ] to form a quadratic equation in terms of the distance variable d v 
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In a first method for obtaining a final solution, the quadratic equation for the 
distance variable d,is solved directly, producing two possible results for the distance 
variable d r This quadratic solution ambiguity can be resolved by discarding a solution 
that is out of the region of interest. The region of interest may be a position fix 
5 determined by the mobile GPS receiver, or may be an indication of the serving cell 

location. In a second method for determining a final solution, the quadratic equation for 
the distance variable d x is solved by processing an additional TOA measurement, i.e., a 
fourth TOA measurement TOA 4 and a third difference TOA 4 -TOA r The fourth TOA 
measurement TOA 4 leads to a third linear equation, similar in form to equations (11) and 
10 (12). This permits a direct, unambiguous solution for x, y and d, (i.e., no approximate 
position knowledge is required). Either final solution method produces an estimated 
range d e to the BS, which may be substituted into equation (6) to compute the BS 
coordinates x, y. 

Tn step 414, the perceived time bias offset b' is determined by subtracting the 

15 estimated range d c (i.e., the calculated distance determined above) from any of the TOA 
equations TOA r In step 416, the computed BS location and the perceived time bias offset 
b' are optimized in order to improve their accuracy and usefulness. 

The nonlinear solution for the BS coordinates x, y provides a starting point for 
calibrating the BS location in the MS. The calculated BS location, however, is capable of 

20 being improved. The calibration may be fine-tuned by linearizing the TOA equations (1)- 
(3) and then computing a Weighted Least Squares (WLS) solution. The WLS solution 
will, in general, improve the accuracy of the calibration of the BS location and will enable 
a reasonable assessment of the accuracy of the calibration. This assessment can be used 
to determine when the BS location calibration is complete (i.e., satisfactorily accurate). 

25 The TOA measurements (TOA,, TOA 2 , TOA 3 ) are linearized about the initial, 

nonlinear, solution by subtracting an estimated TOA (computed from the initial BS 
location with b' set to zero) from each measured TOA. Once linearized about the initial 
solution, the TOA measurement set (TOA 17 TOA,, etc.) can be collected together as a 
vector of measurement residuals, TOA res , defined as: 

30 TOA res = Hx sol (13) 

where H is a measurement gradient vector and x sol is a three dimensional solution vector. 
Each row of H has the following general form: 
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where u x and u y are the components of a Line Of Sight (LOS) unit vector to the BS 
(computed using its initial estimated location). 
The WLS x sol solution is given by: 

Xsol = (H T RWH T R- l TOA rcs (14) 
5 where R is a measurement error covariance matrix corresponding to a TOA 

measurement set. Note that the R matrix must include the effects of the error in the GPS 
solution itself, projected along the estimated LOS unit vector to the BS. 

A by-product of the WLS solution x sol is the solution error covariance' matrix (H T R~ 
l H)'\ which can and should be used to predict the accuracy of the calculated BS location 
10 and the perceived time bias offset b'. 

As previously mentioned, the calibration is complete when the solution error 
variances become acceptably small, and the error variances are confirmed by the 
residuals of an over-determined solution. A proposed acceptable solution error variance 
threshold is one meter, one sigma, which is negligible relative to the positioning 
15 requirements (roughly fifty meters, one sigma, as specified in a U.S. Federal 
Communication Commission (FCC) edict on emergency positioning). As this 
requirement is relaxed, the required calibration accuracy can be relaxed proportionately. 

Upon completion of the BS location calibration, subsequent TOA measurements 
can be combined with GPS pseudorange measurement data to improve the location 
20 accuracy of the MS for calibrating other BS locations. 

If the BS location is already known to the GPS receiver, but the perceived time 
bias offset b' is not, the problem is already in linear form and the solution is 
straightforward. The TOA residual vector is directly related to the time offset bias as: 

TOA=Hb e (15) 

res e 

25 Each row of H is unity and the perceived time bias offset b c collapses to a weighted 

average of the measurement residuals. The accuracy of the perceived time bias offset b c 
can be determined from the solution error covariance matrix (FfR'H)" 1 . 

FIG. 5 shows a flowchart 500 of a second embodiment of a method for calibrating 
BS locations and perceved time bias offsets. In step 501, an MS collects TOA 

30 measurements from a BS in a manner similar to the first flowchart 400, producing 
equations (l)-(3). In step 505, a corresponding number of GPS receiver position 
measurements are generated by the GPS receiver. 



9 
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In step 506, the singular geometry of the TOA measurement set is again tested 
using the matrix M. As before, if the singularity of the received TOA measurements is 
inadequate, the method may acquire more TOA measurements and wait until the 
singularity reaches a certain threshold. 
5 In step 508, the TOA measurements are used to create TOA difference equations 

(4) and (5). In step 511, the difference equations have the distance variable d x added to 
both side of the equation, the equations are squared, and the term d^ is subtracted from 
both sides to yield equations (11) and (12). In step 514, the distance variable d, is 
substituted for the value V[(x - xj 2 + (y - y 4 ) 2 ] to produce a quadratic equation in terms of 

10 the distance variable d r 

In step 520, the distance variable d x is solved. In a first final solution method, the 
quadratic equation for the distance variable d,is solved directly, producing two possible 
results for the distance variable d r This quadratic solution ambiguity can be resolved by 
discarding a solution which is out of the region of interest. The region of interest may be 

15 a position fix determined by the mobile GPS receiver, or may be a set of previously 
acquired BS coordinates. 

In a second final solution method, the quadratic equation for the distance variable 
d, is solved by processing a fourth TOA measurement TOA 4 and a third difference TOA 4 ~ 
TOA r The third difference TOA 4 -TOA 1 leads to a third linear equation, similar in form to 

20 equations (11) and (12). This permits a direct, unambiguous solution for x, y and d, (i.e., 
no approximate position knowledge is required). 

In step 523, the computed BS coordinates x, y may be used to compute the 
perceived time bias offset b'. The perceived time bias offset b' may be determined by 
subtracting the estimated range d e from any of the TOA equations TOA t . In step 527, the 

25 computed BS location and the perceived time bias offset b' are optimized in the same 
manner as in the first flowchart 400. 

Once the BS calibration is essentially complete, the perceived time bias offset b' is 
periodically re-determined and associated with each MS location or BS service area. Over 
time, this association permits construction of a temporal map, or database, of perceived 

30 time bias offsets versus MS position over the area served by a BS. The perceived time 
bias offset may represent the time bias offset b as affected by spatial error from 
multipath in addition to the fixed error induced by cable delays, internal BS harnessing, 
software processing, clock source distribution, and other factors. This map can be 
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derived and maintained for many BSs by the cellular infrastructure and broadcast to 
other MSs (GPS of non-GPS equipped) to improve their location accuracy. Various other 
information can be appended to this map, such as a confidence factor based on a solution 
error variance, a time stamp based on ephemeris data, or the time bias offset b, for any 
5 perceived time bias offset b'. Also the map may be compressed by reducing it merely to 
a correction factor of b. 

The perceived time bias offset may be stored in MSs and subsequent offset 
measurements loaded periodically to a recording device in the infrastructure capable of 
creating the map. Alternatively, a simultaneous transaction session, stich as a short 

10 message service (SMS), Bluetooth, or a defined protocol, could act to transfer the 

perceived time bias offsets over the air back to the infrastructure or another MS, which 
may or may not be GPS assisted. The described perceived time bias offset calculation 
may be realized as a standard mode of operation requested by the infrastructure thereby 
allowing a population of MSs to participate in realizing the perceived time bias offset 

15 map. 

The perceived time bias offset map can be used in a variety of applications to aid 
in more accurate location of an MS, whether or not GPS-equipped. For example, typical 
multipath reduction techniques attempt to remove a calculated time bias offset based on 
a statistical measure of the variance of a signal over time. When line-of-sight (LOS) 

20 conditions exist (i.e., there is little multipath), multipath reduction techniques end up 
compensating for a non-existent problem, and this "compensation" in effect adds error 
to the TOA estimate. To intelligently aid the effectiveness of such multipath reduction 
techniques, the realized map of perceived time bias offset versus position over the area 
served by each BS may be queried. The returned perceived time bias offset can be used 

25 as a discriminator to determine when to use a multipath reduction technique. 

Accessing precise perceived time bias offset map information could be obtained 
once a first GPS-based position is obtained. This precise value acts as a multipath 
reduction technique discriminator for BS TOA signals used in conjunction for 
subsequent hybrid GPS and cellular terrestrial position fixes. Alternatively, after a first 

30 terrestrial-only (enhanced observed time difference (EOTD) or advanced forward-link 
trilateration (AFLT)) location fix is obtained via blindly applying a multipath rejection 
technique, the perceived time bias offset map values corresponding to the sum of the 
delays for each BS could be used in subsequent fixes as a location-specific reduction 
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value for terrestrial-only fixes. This application assumes an MS to also contain terrestrial 
location technology. The effect is to hone the terrestrial measurements over time. 

Approximate perceived time bias offset information could be statistically 
generated for a cell or sector's service area. In rural cells, typical scenarios may incur low 
5 amounts of multipath, in contrast to dense urban environments. This method allows the 
MS to decide if multipath rejection techniques should be executed. The lifetime of the 
specific perceived time bias offset value should ideally be re-evaluated upon each GPS 
position calculation opportunity that results in a BS calibration operation. GPS time, such 
as that found in ephemeris data, can be used to evaluate the "staleness" of the perceived 

10 time bias offset value. 

FIG. 6 shows a GPS receiver 600 that can be used for calibrating BS locations and 
perceived time bias offsets. The GPS receiver 600 includes an antenna 603, a GPS receiver 
component 608, a processor 613, and storage 620. The processor 613 may be any type of 
general purpose processor available in the art. 

15 The storage 620 may be any type of storage memory, including random access 

memory (RAM), magnetic disk or tape, optical memory, etc. The storage 620 includes a 
TOA measurement storage 621 capable of storing a plurality of TOA measurements, a 
position fix storage 623 capable of storing a plurality of position fixes of the GPS receiver 
600, a BS coordinates storage 626 capable of storing a calculated BS location and/or a BS 

20 location received from the infrastructure, a time bias offset storage 628 capable of storing 
a calculated perceived time bias offset b' and/ or a time bias offset b received from the 
infrastructure, an approximate location storage 631 capable of storing an approximate 
geographical location of the BS, a TOA residual vector storage 634 capable of storing a 
TOA residual vector formed from the TOA measurements, and a solution error variance 

25 storage 636 capable of storing a solution error variance threshold. 

It should be noted that although the individual components are shown as being 
dedicated components, alternatively the components may be shared with other 
hardware and functions. For example, when the GPS receiver 600 is combined with a 
cellular handset in an MS, the processor 613 and the storage 620 may be shared between 

30 the cellular phone functionality and the GPS receiver functionality. 

In use, a GPS signal is received through the antenna 603 and decoded by the GPS 
receiver component 608. The GPS receiver component 608 sends the decoded GPS data 
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to the processor 613, which may store some or all of the received GPS data in the storage 
620. 

The processor 613 is capable of collecting and storing TOA measurements, 
collecting and storing a corresponding number of GPS position measurements, 
5 subtracting a first TOA measurement from subsequent TOA measurements to create 

TOA differences, substituting a distance variable d 1 for the unknown BS coordinates x, y 
in the TOA differences to produce a quadratic equation, solving the quadratic equation 
for the distance variable d x to obtain an estimated range d e , solving for the unknown BS 
coordinates x, y, obtaining the perceived time bias offset b' by subtracting the estimated 

10 range d e from the TOA measurements, linearizing TOA equations to produce linearized 
TOA equations, generating a vector of measurement residuals TOA res , computing a 
weighted least squares solution x^ of the TOA equations, using the weighted least 
squares solution x sol to assess an accuracy of the linearized TOA equations, and stopping 
the linearizing when solution error variances of the weighted least squares solution x^ 

15 become less than a solution error variance threshold. This solution error variance 
threshold can either be predetermined or calculated and modified dynamically. 

In this manner, the GPS receiver 600 may operate autonomously and does not 
need to be dependent on reception of a BS location and time bias offset b. The BS location 
and perceived time bias offset b' derived by the GPS receiver 600 may be used in 

20 conjunction with TOA measurements to accurately calculate a position fix of the GPS 
receiver 600. 

While the invention has been described in detail above, the invention is not 
intended to be limited to the specific embodiments as described. For example, satellite 
positioning systems such as GLONASS can be substituted for the GPS satellite 
25 positioning system. It is evident that those skilled in the art may now make numerous 
uses and modifications of and departures from the specific embodiments described 
herein without departing from the inventive concepts. 

What is claimed is: 
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CLAIMS 



1. A method for calibrating a base station location using a satellite positioning 
receiver in a mobile station, characterized by the steps of: 

collecting a number of time-of-arrival (TO A) measurements arriving at the mobile 
station from a base station; 

simultaneously collecting a corresponding number of satellite position 
measurements of the satellite positioning receiver; 

determining that the number of TOA measurements and the corresponding 
number of satellite position measurements is sufficient for calibrating the base station 
location; and 

computing the base station location using the number of TOA measurements and 
the corresponding number of satellite position measurements. 

2. The method of claim 1, further characterized by the step of: 
computing a perceived time bias offset b' using the number of TOA 

measurements and the corresponding number of satellite position measurements. 

15 3. The method of claim 2, further characterized by the steps of: 

associating the perceived time bias offset b' with a satellite positioning receiver 
position; and 

creating at least a portion of a map of perceived time bias offset as a function of 
satellite positioning receiver position. 



20 4. The method of claim 3, further characterized by the step of: 
transmitting the map to another mobile station. 

5. The method of claim 3, wherein the step of creating comprised: 
updating the map within a cellular infrastructure. 

25 6. The method of claim 5, further characterized by the step of: 

using the perceived time bias offset in the map to determine whether to use a 
multipath reduction technique. 
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7. The method of claim 6, further characterized by the step of: 

using the perceived time bias offset in the map to compute a position of the 
mobile station. 

8. The method of claim 1, wherein the step of computing comprises the steps of: 
5 • subtracting a first TOA measurement from subsequent TOA measurements to 

create two TOA differences TOA 31 , TOA 21 ; 

substituting a distance variable dj for base station coordinates x, y in the two TOA 
differences TOA^, TOA^ to produce a quadratic equation for the distance variable d lf 
where the distance variable d, = V[(x-x t ) 2 + (y-y t ) 2 ]; 
10 solving the quadratic equation for the distance variable d x to obtain an estimated 

range d e ; and 

solving for the base station coordinates x, y using the estimated range d e . 

9. The method of claim 8, further characterized by the step of: 
computing a perceived time bias offset b' using the number of TOA 

15 measurements and the corresponding number of satellite position measurements. 

10. The method of claim 9, wherein the step of computing a perceived time bias offset 
b' comprises: 

subtracting the estimated range d e from a TOA measurement. 

11. A method for calibrating a base station location in a mobile station having a 
20 satellite positioning receiver, characterized by the steps of: 

collecting at least three TOA measurements TOA a , TOA 2 , TOA 3 arriving at the 
mobile station from a base station at different times t, wherein a TOA measurement is 
capable of being expressed as TOA t = b' + V[(x - xf + (y - y t ) 2 ], with variable b' being a 
perceived time bias offset of the base station and (x, y) being base station coordinates; 
25 simultaneously collecting at least three satellite position measurements of a 

position of the satellite positioning receiver, wherein a position measurement is denoted 
as x t and y t as a function of time; 
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determining through a geometry matrix that the at least three TOA measurements 
and the at least three satellite position measurements are sufficient for determining the 
base station location; 

subtracting a first TOA measurement TOA x from at least two subsequent TOA 
5 measurements TOA 2/ TOA 3 to create at least two TOA differences TOA 31 , TOA^; 

substituting a distance variable d t for the base station coordinates x, y in the at 
least two TOA differences TOA 31 , TOA^ to produce a quadratic equation for the distance 
variable d u where the distance variable d, = V[(x-x t ) 2 + (y-yf]; 

solving the quadratic equation for the distance variable d t to obtain an estimated 
10 range d e ; 

solving for the base station coordinates x, y using the estimated range d e ; and 
obtaining the perceived time bias offset b' by subtracting the estimated range d e 
from a TOA measurement. 

12. The method of claim 11, wherein the step of solving the quadratic equation 
15 comprises the steps of: 

solving the quadratic equation to produce two solutions; and 
discarding a solution that is out of a geographical region of interest to leave an 
accepted solution. 

13. The method of claim 11, wherein the step of solving the quadratic equation 
20 comprises the steps of: 

obtaining at least a fourth TOA measurement and at least a third TOA difference; 

and 

using the at least a third TOA difference in conjunction with the quadratic 
equation to determine a solution for the distance variable d x and for the base station 
25 coordinates x, y. 

14. The method of claim 11, further including the steps of: 

linearizing the at least three TOA measurements TOA u TOA 2 , TOA 3 by 
subtracting estimated ranges from them, with the estimated ranges determined from the 
base station coordinates x, y, to produce linearized TOA equations; 
30 generating a vector of measurement residuals TOA res ; 
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computing a weighted least squares solution of the linearized TOA equations; 
correcting the base station coordinates x, y using the weighted least squares 
solution x sol ; 

repeating the steps of generating, computing, and correcting until a magnitude of 
5 the weighted least squares solution x^ becomes less than a solution error variance 
threshold; 

using a covariance matrix of the weighted least squares solution x sol to assess an 
accuracy of the base station coordinates x, y; and 

using the accuracy to determine that the base station location calibration is 
10 complete. 

15. The method of claim 14, wherein the step of using the covariance matrix 
additionally makes use of measurement residual information from the weighted least 
squares solution x sol . 

16. The method of claim 14, wherein the step of using the accuracy additionally 
15 makes use of residual information from the weighted least squares solution x^. 

17. The method of claim 14, wherein the vector of measurement residuals TOA res is 
given by TOA res = H x sol , where H is a measurement gradient vector and the weighted 
least squares solution x sol is a three dimensional weighted least squares solution vector, 
and where a row of the measurement gradient vector H has a form of h* = [ u x u y 1] with 

20 u x . and u y being components of a line of sight unit vector to the base station. 

18. The method of claim 14, wherein the weighted least squares solution x sol is given 
as x SOI = (H^HyWR^TOA^, where R is the covariance matrix. 

19. The method of claim 14, wherein if the base station coordinates x, y are known but 
the perceived time bias offset b' is not known, then a TOA residual vector TOA res is given 

25 by TOA res = H b e , where each row of H is unity and an estimated time bias offset b e 
collapses to a weighted average of measurement residuals, and an accuracy of the 
estimated time bias offset b e can be determined from the covariance matrix. 
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20. An apparatus for calibrating a base station location in a mobile station, 
characterized by: 

a satellite positioning receiver; 
a radio transceiver; 

5 a time-of-arrival (TO A) measurement storage for storing a plurality of TOA 

measurements; 

a satellite positioning receiver position fix storage for storing a plurality of 
position fixes; 

a base station coordinates storage; 

10 a time bias offset storage; and 

a processor, for communicating with the radio transceiver and satellite 
positioning receiver to receive signals from satellite positioning satellites and from a base 
station, and further communicating with the TOA measurement storage, the satellite 
positioning receiver position fix storage, the base station coordinates storage, and the 

15 time bias offset storage, with the processor capable of collecting and storing TOA 
measurements, collecting and storing a corresponding number of satellite position 
measurements, subtracting a first TOA measurement from subsequent TOA 
measurements to create TOA differences, substituting a distance variable d, for base 
station coordinates x, y in the TOA differences to produce a quadratic equation, solving 

20 the quadratic equation for the distance variable d a to obtain an estimated range d e , 

solving for the base station coordinates x, y, obtaining a perceived time bias offset b' by 
subtracting an estimated range d e from the TOA measurements, linearizing the TOA 
measurements to produce linearized TOA equations, generating a vector of 
measurement residuals TOA res , computing a weighted least squares solution x sol of the 

25 linearized TOA equations, using the weighted least squares solution x^ to assess an 

accuracy of the linearized TOA equations, and repeating the generating, computing, and 
using until a solution error variance of the weighted least squares solution x sol becomes 
less than a solution error variance threshold. 

21. The apparatus of claim 20, further including an approximate geographical 
30 location storage, communicating with the processor, for storing an approximate 

geographical location for use in obtaining an estimated range d e . 
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22. The apparatus of claim 20, further including a TOA residual vector storage, 
communicating with the processor, for holding the vector of measurement residuals 
TOA rcs . 

23. The apparatus of claim 20, further including a solution error variance threshold 
5 storage, communicating with the processor, for holding the solution error variance 

threshold. 

24. A map characterized by: 

a perceived time bias offset as seen by a mobile station as a function of a 
computed satellite positioning receiver position. 

10 25. A map according to claim 24 further characterized by: 

a confidence factor associated with the perceived time bias offset. 

26. A map according to claim 24 further characterized by: 

a time stamp associated with the perceived time bias offset. 

27. A map according to claim 24 further characterized by: 

15 a time bias offset associated with the perceived time bias offset. 
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